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ABSTRACT 

Aims. This paper aims to study the configuration of two components caused by rotational and tidal distortions in the model of a binary 
system. 

Methods. The potentials of the two distorted components can be approximated to 2nd-degree harmonics. Furthermore, both the 
accretion luminosity (cr,) and the irradiative luminosity are included in stellar structure equations. 

Results. The equilibrium structure of rotationally and tidally distorted star is exactly a triaxial ellipsoids. A formula describing the 
isobars is presented, and the rotational velocity and the gravitational acceleration at the primary surface simulated. The results show 
the distortion at the outer layers of the primary increases with temporal variation and system evolution. Besides, it was observed 
that the luminosity accretion is unstable, and the curve of the energy-generation rate fluctuates after the main sequence in rotation 
sequences. The luminosity in rotation sequences is slightly weaker than that in non-rotation sequences. As a result, the volume expands 
slowly. Polar ejection is intensified by the tidal effect. The ejection of an equatorial ring may be favoured by both the opacity effect 
and the g e (9, <^)-effect in the binary system. 
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1. Introduction 

In the conventional model of binary stars, there is no con- 
sideration of spin and tidal effects (Eggleton 1971,1972,1973; 
Hofmeister, Kippenhahn & Weigert, 1964; Kippenhahn et al. 
1967; etc.); however, rotation and tide have been regarded as 
two important physical factors in recent years, so they need to be 
considered for a better understanding of the evolution of massive 
close binaries (e.g., Heger, Langer & Woosley 2000a; Meynet 
& Maeder 2000). The structure and evolution of rotating sin- 
gle stars has been studied by many investigators (Kippenhahn 
& Thomas [19701 E ndal & Sofia [19761 Pinsonneaul et al. [19891 
Meynet & Maeder [19971 Langer 1998, 1999; Huang 2004a). 
However, it is also very important to study the evolution of ro- 
tating binary stars (Jackson 1970; Chan & Chau 1979; Langer 
2003; Huang 2004b; Petrovic et al 2005a,b; Yoon et al. 2006). 
The effect of spin on structure equations has been investigated 
(e.g. the present Eggleton's stellar evolution code; Li et al 
2004a,b, 2005; Kahler 2002). They adopted the lowest-order 
approximate analysis in which two components were treated as 
spherical stars. In fact, with the joint effects of spin and tide, the 
structure of a star changes from spherically symmetric to non- 
spherically symmetric. Then, the stellar structure equations be- 
come three dimensional. Theory distinguishes two components 
in the tide, namely equilibrium tide (Zahn 1966) and the dynam- 
ical tide (Zahn 1975). Then, the dissipation mechanisms acting 
on those tides, namely the viscous friction for the equilibrium 
tide and the radiative damping for the dynamical tide, have been 
identified (Zahn 1966, 1975, 1977). The distortion throughout 
the outer regions of the two components is not small in short- 
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period binary systems. The higher-order terms in the external 
gravitational field should not be ignored (Jackson 1970). 

It is a very complex process to determine the equilibrium 
structure of the two components. Therefore, approximate meth- 
ods have been widely adopted for studying these effects. In 
1933, the theory of distorted polytropes was introduced by 
Chandrasekhar. Kopal (1972,1974) developed the concept of 
Roche equipotential and of Roche coordinates to analyse the 
problem of rotationally and tidally distorted stars in a binary 
system. Bursa (1989a,1988) took advantage of the high-order 
perturbing potential to describe rotational and tidal deforma- 
tions to discuss the figures and dynamic parameters of syn- 
chronously orbiting satellites in the solar system. The equilib- 
rium structure of the two components were treated as two non- 
symmetric rotational ellipsoids with two different semi-major 
axes a\ and (a\ > ai) by Huang (2004b). It is very important 
that Kippenhahn & Thomas (1970) introduced a method of sim- 
plifying the two-dimensional model with conservative rotation 
and allowed the structure equations for a one-dimensional star 
to incorporate the hydrostatic effect of rotation. This method has 
been adopted by Endal & Sofia (1976) and Meynet & Maeder 
(1997), who applied it to the case of shellular rotation law (Zahn 
1992). In this case, the rotation rate takes the simplified form of 
Q. = Q(r). It was demonstrated that the shape of an isobar in the 
case of the shellular rotation law is identical to one of the equipo- 
tentials in the conservative case of Meynet & Maeder (1997). 

At the semi-detached stage, both mass transfer between the 
components and luminosity change of a secondary exist due to 
the release of accretion energy which is correlative with the ex- 
ternal potential of the two components. When the joint effect of 
rotation and tide are considered, the potential of the two compo- 
nents are different from those in non-rotational cases. Therefore, 
the luminosity due to the release of accretion energy, as well 
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as irradiation energy, can significantly alter the structure and 
evolution of the secondary. In a rotating star, meridional circu- 
lation and shear turbulence exist, both of which can drive the 
transport of chemical elements. This effect is stronger and has 
already been studied by many scholars (Endal & Sofia 1978; 
Pinsonneaul et al. [T9891 Chaboyer & Zahn 1992; Zahn 1992; 
Meynet & Maeder 1997; Maeder 1997; Meader & Zahn 1998; 
Maeder & Meynet 2000; Denissenkov et al. 1999; Talon et al. 
1997; Decressin et al.2009). In this paper, the amplitude expres- 
sion for the radial component of the meridional circulation ve- 
locity U (r) considers the effect of tidal force, which may be im- 
portant in a massive close binary system. 

This paper is divided into four main sections. In section 2, 
the structure equations of rotating binary stars are presented. 
Material diffusion equations and boundary conditions are pro- 
vided. Then, the accretion luminosity, including gravitational en- 
ergy, heat energy, and radiation energy, is deduced. In section 3, 
the results of numerical calculation are described and discussed 
in detail. In section 4, conclusions are drawn. 

2. Model for rotating binary stars 

2.1. Potential of rotating binary stars 

It is well known that the rotation of a component is synchronous 
with the orbital motion of a system thanks to a strong tidal ef- 
fect. Such synchronous rotation also exists inside the component 
(Giuricin et al. I1984I Van Hamme & Wilson 11990b ; therefore, 
conventional theories usually assume that two components ro- 
tate synchronously and revolve in circular orbits (Kippenhahn & 
Weige rt [T967l De Loo re [T980l Huang & Taam [19Ml Vanbeveren 
1991 ; De Greve 1993 ). A coordinate system rotating with the or- 
bital angular velocity of the stars is introduced. The mass centre 
of the primary is regarded as the origin, and it is presumed that 
the z-axis is perpendicular to the orbital plane, and the positive 
x-axis penetrates the mass centre of the secondary. The gravita- 
tional potential at any point P(r, 9, if) of the surface of the pri- 
mary can be approximately expressed as 



¥ = V + -OVfl 



P 2 (cos0)) + V„ 



(1) 



where V is the gravitational potential and given by Bur.se 
(1989a,1988), 



V = 
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Here, V, is the tidal potential (Burse 1989a) 
GM 2 
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(3) 



where it is assumed that the mean equatorial radius equals that of 
the equivalent sphere in the above equation for the convenience 
of calculation. Both M\ and M 2 are the mass of the primary and 
the secondary, respectively, and r p represents each equivalent ra- 
dius inside the star, P^cosO) and P^icosO) are the associated 
Legendre function(P°(cos 0) = § cos 2 6-\, P 2 (cos 0) = 3 sin 2 6), 
D is the distance between the two components, and Q. is the or- 
bital angular velocity of the system. It can be represented by 

O 2 = G(Mi + M 2 )/D 3 , (4) 

where 7* 0) and are dimensionless stokes parameters. If M\ 
can generally be negligible compared to M 2 , the stokes parame- 
ters can be expressed as (Burse 1989a,1988) 
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where k s is the secular Love number, which is expressed as a 
measure of the body-yield-to-centrifugal deformation, and k, is 
an analogous parameter that is introduced to describe the secu- 
lar tidal deformations. The response of the body to its centrifu- 
gal acceleration and to the tidal perturbing potential is different 
in the usual case. Therefore, the body-yield-to centrifugal de- 
formation is not equal to the body-yield-to-tidal deformation. 
If the subject investigated is regarded as an ideal elastic body, 
the body-yield-to centrifugal deformation is equal to the body- 
yield-to-tidal deformation, k s - k,. In the ideal static equilib- 
rium, k s = k, — 1 (Burse 1989a). We assume the ideal static 
equilibrium in this paper, q is the mass ratio of the secondary 
to the primary (q = ^). With Eqs. (2) and (3) being combined 
with Eq. (1), the potential of the primary can be obtained as 



GM\ D 1 D (r p \ 
= -=!{- + -- - [-7 2 (3cos 2 0- 1) 



D r\ 2 r\\r\ 

,<2) 



(7) 
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+6J\ ' sin z 6» cos 2^] + -(1 + ?X^) sirr 
+ i^) 2 [3sin 2 0(l+ CO s2^)-2]}, 

The potential of the secondary is deduced by substituting M 2 
for Mi and - for q. The isobar defined by the equation P - 
const is assumed to be a triaxial ellipsoid with three semi-major 
axes: a, b, and c. The shortest axis defined by c is identical to its 
rotational axis and perpendicular to its orbital plane. The longest 
axis defined by a is identical with its x-axis. 



2.2. Considering stellar structure equations with spin and 
tidal effects 

The spin of the two components is rigid rotation, and it belongs 
to conservative rotation. The definition of equivalent sphere was 
adopted in a practical calculation. Therefore, the triaxial ellip- 
soid model is simplified to a one-dimensional model. The struc- 
ture equations are presented as 
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where <x flc is the energy source per unit mass caused by mass 
overflow and irradiation. Because accretion luminosity is caused 
by energy sources in the gainer's outermost layer, there exists 



cr nr Am = AL a , 



(11) 



where Am is the photosphere mass of the secondary. The sur- 
face temperature of the secondary may be approximated by the 
formula, L 2 + AL acc = AnBfyrT\p where L 2 is the luminosity 
coming to the photosphere from the stellar interior, and <x is the 
Stefan-Boltzmann constant: 



dlnT 
(5) d\nP 
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(15) 



where < g e ff > and < g~j^ > are the mean values of effective 
gravity and its opposites over the isobar surface, and V# is the 
radiative temperature gradient. The factors f P and f T depend on 
the shape of the isobars. 

2.3. Calculation of quantities f P and f T 

2.3.1 . Shape and gravitational acceleration of triaxial 
ellipsoid 

To obtain the factors f P and f T , the mean values < g e ff > and 
< g~jf > over the isobar surface have to be calculated. Therefore, 
the shape of isobars must be given first. The functions for the 
semi-major axes a, b, and c to the radius of the equivalent sphere 
r P can be obtained from Eq. (7) as 
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The left hand side of Eq.(17) corresponds to 6 = | and tp = 0, 
while the one of Eq.(18) corresponds to 6 = | and tp = | . The 
three semi-major axes a, />, and c of a triaxial ellipsoid can be 
obtained numerically by solving (16), (17), and (18). From (7), 
the quantities g r , g e , and g v at the surface of the two components 
take the forms of 



gr dr 
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the cylindrically symmetric potential due to rotation; \p t the 
non-symmetric potential due to tidal force, and fa, the non- 
symmetric part of the gravitational potential due to the distor- 
tion of the component considering the rotational and tidal effects. 
Therefore, the total potential at P(r, 9, ip) is 

m = fa+fa+^+fa 
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The quantity 77, can be evaluated by numerically integrating the 
Radau's equation (cf. Kopal 1959) 



dVj pirn) 
ro-± + 6^(1/; ; + l) + 77/77, 
</r p(r ) 



1) = 7(7 + 1). 
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for j=2,3,4, and boundary condition 77/0) = j - 2. The quan- 
tity ro is the mean radius of the corresponding isobar. The local 
effective gravity is given by differentiation of the total potential 
and is written as 



, 1 = [ (^) 2 + ( 1 ^ 2 



r + (- 



1 = 1,2). 



(24) 



dr r 80 r sin 6 dtp 

The integral in above equations and their derivatives must be 
evaluated numerically. The mean values of g e ff t and g~jj, f 
the surfaces of the triaxial ellipsoids can be obtained as 

1 r*^R 

<?effi>=— I I gefftrf sin 6d0d<p,(i = 1,2) 
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According to Eqs. (13) and (14), the values of f P and /j- can 
be obtained when the mean values < g e ff > and < g~^ > are 
known. r\ and r 2 are the distances between the centre of the com- 
ponents and the surfaces of two triaxial ellipsoids. They are 



r 2 = 



a]b\c\ 



b 2 c 2 sin 6 cos 2 tp + a 2 c 2 sin #sin <p + a 2 b 2 cos 2 



The surface area of the isobar can be expressed as 

S p = y(fl 2 + b 2 + C 2 ). 



(27) 
1,2). 

(28) 
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(20) 2.4. Element diffusion process 



(21) 



However, the total potential in the stellar interior (to first- 
order approximation) can be composed by four parts (Kopal 
1959, 1960, 1974; Endal & Sofia 1976 and Landin 2009): fa, 
the spherical symmetric part of the gravitational potential; fa, 



The effect of meridian circulation can drive the transport of 
chemical elements and angular momentum in rotating stars. For 
the components in solid-body rotation, no differential rotation 
exists that can cause shear turbulence. According to Endal & 
Sofia (1978) and Pinsonneault (1989), the transport of chemical 
composition is treated as a diffusion process. The equation takes 
the form of (Chaboyer & Zahn 1992) 



dy a 
dt 



J_d_ 

pr 2 dr 
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where (^-)„« c is a source term from nuclear reactions, and y a 
is the relative abundance of a - th nuclide. The diffusion coef- 
ficient Ddif given by Heger, Langer & Woosley (2000a) can be 
expressed as 



D dif = min{di„ s ,,H v ,EsW(r), 



(30) 



where di„ st and H v es denote the extent of the instability and the 
velocity scale height, respectively. The expression for the ampli- 
tude of the radial component of the meridional circulation ve- 
locity U(r) (derived from Kippenhahn & Weight 1990) has been 
modified to take the effects of radiation pressure and tidal force 
into-account, which are important in a massive close binary sys- 
tem. It is noticed that 



U(r) = 

U V> 3 ^ GM r 



GM 2 r 



L r y-l 1 



GM r D 3 ' grM r y V orf -V 



(31) 



(1 _ Jt -£-) 



where the term + -jj^r is the local ratio of centrifugal force 

and tidal force to gravity, y is the ratio of the specific heats 
L r represents the luminosity at radius r, M r is the mass enclosed 
within a sphere of radius r, V is the actual gradient and V fl ^ is 
adiabatic temperature gradient, s m = — gives the mean energy 
production rate, and e is local generation rate of nuclear-energy. 

There is no source or sink at the inner and the outer bound- 
aries of the two components. Therefore, the boundary conditions 
are used as 



dy a 

(-7— )i=\ 
or 







dy a 
or 



(32) 



where the subscript ; denotes different layers inside stars. The 
initial abundance equals the one at the zero-age main sequence. 
Therefore, the initial condition is 



(y a )i 



i\t=0 



(y a )i 



(33) 



2.5. Luminosity accretion 



In the case where the joint effect of rotation and tide is ignored, 
the two components are spherically symmetric. The star fills 
its Roche lobe and begins to transfer matter to the companion. 
However, in the case with the effects of rotation and tide being 
considered, the components are triaxial ellipsoids. The condition 
for the mass overflow through Roche lobe flow should be revised 
as oi = r rohe (Huang 2004b). It is assumed that the transferred 
mass is distributed within a thin shell at the surface of the pri- 
mary before the transfer, and within a thin shell at the surface 
of the secondary after the transfer. Three forms of energy (in- 
cluding potential energy, heat energy, and radiative energy) are 
transferred to the secondary. The mass transfer rate is m. Two 
different cases are considered: 

a) If the joint effect of rotation and tide is ignored, the accre- 
tion luminosity can be expressed directly in terms of the Roche 
lobe potential at the inner Lagrangian point, } ¥ Ll , and at the sur- 
face of the secondary Y s (Han & Webbink 1999): 



2 ^ D 



g_\2 _ q 

l+q> R 2 /D 



1 



X L JD l-X Ll /D 
l+q. 



1-R 2 /D 



2 Vi D l+q> 



(34) 



where X L \ is the distance between the primary and L\, and R2 is 
the radius of the secondary. 

b) If the joint effect of rotation and tide is considered, the 
equilibrium structure of the two components will be treated as 



triaxial ellipsoids. The release of potential energy because of the 
accretion of a mass rate m to the secondary is given by 

AL P = MV Ll - V S ) = *PIt£jb + T=ijo + 

- ^[f + ||(|) 2 (2/ 2 °) - \q{%)\ (35) 

where *F. S is the potential of the secondary. Similarly, as the two 
components have different temperatures, the transmitted thermal 
energy will be 



3kT, 



effl 



3kT, 



eff2 



2/j.irrip 2fi 2 m p 



). 



(36) 



where T e ff \ and T e ffi represent the effective temperature of the 
primary and the secondary, respectively, and yi\ and //2 are the 
mean molecular weights of the primary and the secondary, re- 
spectively. m p refers to proton mass. Because of the irradiation, 
energy accumulated by the primary and the secondary can be 
given by (Huang & Taam 1990) 



1 



AL r , u = -[1-(D Z -Ri 2 )i/D]l4,u 



(37) 



where Ri and R2 are the radii of the primary and the secondary, 
and L\ and L2 are the luminosities of the primary and the sec- 
ondary, respectively. The total accretion luminosity is 



AL acc = f3AL tot = 0(AL P + AL T + AL r ). 



(38) 



Because a part of the total energy may be dissipated dynamically, 
P is assumed to range from 0.1 to 0.5 (Huang 1993). A value 
P = 0.3 is adopted. 

3. Results of numerical calculation 

The structure and evolution of binary system was traced with 
the modified version of a stellar structure program, which was 
developed by Kippenhahn et al., (1967) and has been updated 
to include mass and energy transfer processes. The calculation 
method is based on the technique of Kippenhahn and Thomas 
(1970) and takes advantage of the concept of isobar (Zahn 1992, 
Meynet and Maeder 1997). Both components of the binary are 
calculated simultaneously. The initial mass of the system com- 
ponents is set at 9M Q and 6M Q . The initial chemical compo- 
sition X equals X=0.70, and Z=0.02 is adopted for the two 
components. Similarly, the initial orbital separation between the 
two components for all sequences is defined as 20.77 1R Q , so 
mass transfer via Roche lobe occurs in case A (at the central 
hydrogen-burning phase of the primary). Two evolutionary se- 
quences corresponding to the evolution with the joint effect of 
rotation and tide being considered or ignored are calculated. The 
sequence denoted by case 1 represents the evolution without 
the effects of rotation and tide being considered, while the se- 
quence denoted by case 2 represents the evolution with the ef- 
fects of rotation and tide being considered. The calculation of 
Roche lobe is taken from the study by Huang & Taam (1990). 
The non-conservative evolution in the two cases was considered. 
Because the local flux at colatitude 9 is proportional to the effec- 
tive gravity g e according to Von Zeipel theorem (Maeder 1999), 
the mass-loss rate due to the stellar winds intensified by tidal, 
rotational, and irradiative effects is obtained according to Huang 
& Taam(1990) (cf. Table 1). The angular velocity of the system 
and the orbital separation between the two components change 
due to a number of factors: changes in physical processes as 
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the binary system evolves, including the loss of mass and an- 
gular momentum via stellar winds, mass transfer via Roche lobe 
overflow, exchange of angular momentum between component 
rotation and the orbital motion of the system caused by tidal ef- 
fect, and changes in moments of inertia of the components. The 
changes in the angular velocity of the system and the orbital sep- 
aration between the two components can be calculated according 
to Huang & Taam (1990), and the results are listed in Table 1. 
Other parameters are treated in the same way for two sequences. 

The evolution of the binary system proceeded as follows (cf. 
Table 1). Evolutionary time, orbital period, mass of two stars, 
luminosities and effective temperature of two stars, central and 
surface helium mass fraction of the primary, and mean equato- 
rial rotational velocities of two stars are listed in Table 1 . Points 
a, b, c, d, e, and f denote the zero-age main sequence, the begin- 
ning of the mass transfer stage, the beginning of H-shell burn- 
ing, the end of central hydrogen-burning, the beginning of the 
central helium-burning stage, and the end of calculation, respec- 
tively. At the beginning of mass exchange, the luminosity and ef- 
fective temperature of the primary component decrease rapidly. 
The secondary accretes 6.174M for case 1 and 5.502M o for 
case 2 during the mass transfer in case A. Because of this mass 
gain, the luminosity and the temperature of the secondary go up. 
When the mass is transferred from the more massive star to the 
less massive one, the separation between the centres of the two 
components as well as the orbital period of the system decrease. 
Some orbital angular momentum is transformed into the spin an- 
gular momentum of both components, and this process is crucial 
to model the spin-up of the accretion star. With mass overflow, 
the mass of the primary will be less than that of the secondary. 
When the mass is transferred from the less massive star to the 
more massive one, the separation between the centres of the two 
components as well as the orbital period of the system increases. 
Some spin angular momenta in both of the components are trans- 
formed into orbital angular momentum. This physical process 
results in a longer epilogue after mass transfer. 

The equilibrium configuration deviates from spherical sym- 
metry because of the centrifugal forces and tidal forces. And the 
deviated region mainly lies in the outer layer of a star. In fact, the 
distorted stellar surface forms the shape of a triaxial ellipsoid. A 
distorted isobar surface can be expressed as 

r = r p [l + f(r)P 2 (cos 0) + g(r)P 2 (cos 6)cos2<p], (39) 

which corresponds to the form of the disturbing potential (Zahn 
1992). The coefficients f(r) and g(r) can be defined as f(r) = 

-% ~ |f and 8(r) = HI - The quantity q is the mean den- 
sity of a star with the mass of M\. It was noticed that at the 
central hydrogen-burning phase, two parameters C\ and C2 in 
Eq. (39) gain the values of 0.703 + 0.125 and 0.491 + 0.102, 
respectively. This formula indicates that the shapes of the two 
components vary with the potentials of the centrifugal force and 
the tidal force. The radial deformation is inversely proportional 
to the mean density of the component. In order to describe the 
distortion, the distribution of the surface rotating velocities of 
the primary is illustrated in Fig. 1 . The four panels (a), (b), (c), 
and (d) correspond to the evolutive time of 0, 2.3386 x 10 7 , 
2.6194 x 10 7 , 2.6287 x 10 7 years and corresponding periods of 
2.776, 2.760, 2.746, and 2.628 days, respectively. The rotational 
velocity rates for the peaks of the semi-major axes b and a are 
t = ZL = 0.9867,0.9401,0.8814, and 0.8664 in four panels, 
respectively. The results show that the surface deformation is 
intensified with the evolution and volume-expansion of the pri- 
mary. The distortion throughout the outer region of the primary 



is considerable. The detailed theoretical models that focus on in- 
vestigation of the outer regions have somewhat deviated from the 
Roche model. The high-order perturbed potential is required for 
studying the structure and evolution of short-period binary sys- 
tems. Matthews & Mathieu (1992) examined 62 spectroscopic 
binaries with A-type primaries and orbital periods less than 100 
days. They concluded that all systems with orbital periods less 
than or equal to three days have circular orbits or nearly circu- 
lar orbits. Zahn (1977) and Rieutord & Zahn (1997) have shown 
how binary synchronization and circularization result from tidal 
dissipation. Based on smoothed particle hydrodynamics (SPH) 
simulation, Renvoize et al.(2002) have quantified the geometri- 
cal distortion effect due to the tidal and rotational forces acted on 
the polytropic secondaries of semi-detached binaries. They sug- 
gest that the tidal and rotational distortion on the secondary may 
not be negligible, for it may reach observable levels of ~ 10% 
on the radius in specific cases of polytropic index and mass ratio. 
Georgy et al.(2008) display that various effects of the rotation on 
the surface of a 2OM star at a metallicity of 10~ 5 and at ~ 95% 
of the critical rotation velocity. They point out that the star be- 
comes oblate with an equatorial- to-polar radius ratio ^ 1.3. 
These results agree closely with ours. 

The variation relative gravitational accelerations, the tidal 
force, and the ratio of f cen /ftid on the surface of the primary un- 
der the coordinate 6 and ip at the beginning of mass overflow 
are shown in Fig. 2. The quantities g r , gg, and g^ are the three 
components of gravitational acceleration. The six panels (a), (b), 
(c), (d), (e), and (f) represent the distribution of g r /g, gelg, giplg, 
gtot/g, ftid, and fcen/ftid, respectively. The quantity g equals the 
gravitational acceleration of the corresponding equivalent sphere 

(g = ^r-)- When the joint effect of rotation and tide is consid- 

r p 

ered, the gravitational accelerations are different from those in 
the conventional model. Gravitational acceleration generally has 
three components. 

It is shown in panel (a) that the relative quantity j reaches 
the maximum value of 1 .048 at the two polar points and drops 
to the minimum value 0.6987 on the equatorial plane because 
the inward tidal force acts on the primary and causes the po- 
lar radius to become shorter. The tidal and centrifugal forces 
pull the primary outwards and change gravitational accelera- 
tions greatly on the equatorial plane. Furthermore, the maxi- 
mum value is 0.9486 and the minimum value is 0.6987 on the 
equatorial plane. The lower values are at the peak of the longest 
axis a and the higher values are at the peak of the axis b. The 
relative quantity j reaches the maximum value of 0.20399 at 

the point of 6 = y + | ; <p = kn, k = 0,1 and vanishes at 
the two polar points and on the equatorial plane in panel (b). 
It can be seen that a secondary maximal value of 0.10214 ex- 
ists at point of = ^ + f; <p = kn + §, k = 0, 1. The rela- 
tive quantity y reaches the maximum value of 0.1050 at point 
= \;<p = y + j, k=0,l,2,3 and decreases to zero at the 
point of <p = y< k = 0, 1,2,3 in panel (c). The total gravi- 
tational acceleration at the surface of the primary is shown in 
panel (d). Its distribution is similar to the one of j because the 
radial component is the maximum value. It is noticed that, as 
expected, the average gravitational acceleration of the rotating 
model is less than for the non-rotating model. It can be observed 
that the tidal force reaches the highest value of 340.05cm/s 2 
at the point of 8 = |;<£ = kn, k = 0, 1 and decreases to the 
lowest of 136.24cm/s 2 at the two polar points. The quantity 
fcenl ftid reaches the maximum value of 2.4905 at the point of 
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Table 1. Parameters at different evolutionary points a, b, c, d, e, and f in sequences of cases 1 and 2. 



Sequence 


Time 
10V 


P 

day 


Mi 
M 


M 2 
M Q 


logLi/La 


logT^ff 


logL 2 /Ls 


logT 2 ,eff 






V rat ,l 
km 1 sec 


V rat ,2 
km 1 sec 


a 

Case 1 
Case 2 


0.0000 
0.0000 


2.777 
2.776 


9.000 
9.000 


6.000 
6.000 


3.639 
3.629 


4.385 
4.381 


3.077 
3.055 


4.287 
4.280 


0.2800 
0.2800 


0.2800 
0.2800 


68.66 


56.39 


b 

Case 1 
Case 2 


2.6725 
2.6267 


2.737 
2.743 


8.929 
8.939 


5.994 
5.995 


3.959 
3.914 


4.285 
4.295 


3.110 
3.085 


4.267 
4.265 


0.8744 
0.8247 


0.2800 
0.2800 


143.69 


63.39 


c 

Case 1 
Case 2 


2.6854 
2.6329 


3.190 
3.192 


5.314 
5.318 


9.608 
9.616 


3.457 
3.264 


4.120 
4.135 


3.811 
3.834 


4.396 
4.406 


0.8799 
0.8267 


0.2801 
0.2801 


121.81 


67.59 


Cl 

Case 1 
Case 2 


3.0784 
3.3567 


11.743 
7.213 


2.727 
3.388 


12.168 
11.497 


3.776 
3.442 


4.128 
4.162 


4.145 
4.123 


4.461 
4.405 


0.9800 
0.9800 


0.5892 
0.3481 


58.41 


41.79 


e 

Case 1 
Case 2 


3.0941 


32.534 


1.797 


13.082 


3.959 


4.028 


4.276 


4.539 


0.9799 


0.8775 






f 

Case 1 
Case 2 


3.1484 
3.5470 


42.382 
32.531 


1.572 
1.714 


13.245 
13.037 


3.186 
3.280 


4.775 
4.686 


4.308 
4.308 


4.549 
4.419 


0.8856 
0.9800 


0.8793 
0.8261 


0.97 


10.74 




= |; (p = kn + |, k = 0, 1 and reaches the secondary maximal 
value of 1.2445 at the point of 6 = f; <p = kn, k = 0, 1. The 
results show that the effect produced by tidal distortion is lower 
in comparison with what is produced by rotational distortion on 
the equatorial plane. However, with the mass conversion, the op- 



posite situation can emerge. It is concluded that tidal distortions 
are related to the mass ratio of the secondary to the primary. 
These results suggest that rotation and tide have strong influ- 
ences on the stellar surface. They modify the gravity and change 
the spherically-symmetric shape into the triaxial ellipsoid shape. 
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Fig. 2. Variation of relative gravitational accelerations at the surface of primary under coordinate 6 and tp as mass overflow begins. 
The quantities g n g e , and g^ are the three components of the gravitational acceleration g tot . Quantity g equals the gravitational 
acceleration of the corresponding equivalent sphere (g = ^r 1 ). 

fry 



Furthermore, the stellar structure equations are basically revised 
due to the distribution of the relative quantity in the outer region. 

According to the Von Zeipel theorem, the mass loss due to 
stellar winds should be proportional to local effective gravity. 
Polar ejection is intensified by the tidal effect. The higher grav- 
ity at the peak of the axis b makes it hotter. The ejection of an 
equatorial ring may be favoured by both the opacity effect and 
the higher temperature at the peak of the semi-axis b. This ef- 
fect is called the g e {6, y?)-effect in this paper. It is predicted that 
the g e (9, ^-effect is as important as the g e -effect suggested by 
Maeder (1999) and Maeder & Desjacques (2001). The shapes 
of planetary nebulae that deviate from spherical symmetry (ax- 



isymmetrical one in particular) are often ascribed to rotation or 
tidal interaction (Soker 1997). Frankowski and Tylenda (2001) 
suggest that a mass-losing star can be noticeably distorted by 
tidal forces, thus the wind will exhibit an intrinsic directivity 
and may be globally intensified. Interestingly enough, the group 
of the B [e] stars shows a two-component stellar wind with a hot, 
highly ionized, fast wind at the poles and a slow, dense, disk-like 
wind at the equator (Zickgraf 1999). Maeder and Desjacques 
(2001) have noticed that the polar lobes and skirt in rj Carinae 
and other LBV stars may naturally result from the g e ff and k- 
effects. Langeret al. (1999) have shown that giant LB V outbursts 
depend on the initial rotationrate. Tout and Eggleton (1988) pro- 
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Fig. 3. Time variation of relative accretion luminosity at semi-detached stage. Panel (a) represents case 1 and panel (b) represents 
case 2. The solid, dotted, dashed and dotted-dashed curves correspond to the relative accretion luminosity with respect to total, 
thermal, potential and irradiative energies, respectively. 
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Fig. 4. Panel (a): variation in total H-burning generation energy rate in two cases. Panel (b): time variation in total H-burning 
generation energy rate in case 2 after main sequence. The solid curve represents case 2 and the dashed curve represents case 1 . 



posed a formula according to which the tidal torque would en- 
hance the mass-loss rate by a factor of 1 + B x (^) 6 , where 
B is a parameter free to be adjusted (ranging from 5 x 10 2 to 
10 4 ). Mass loss and associated loss of angular momentum are 
anisotropic in rotating binary stars. The theories for describing 
the mass loss and angular momentum loss from stellar winds 
should be altered partly in future work. 

The time variation of relative accretion luminosity at the 
semi-detached stage is shown in Fig. 3. The two panels (a) and 
(b) correspond to cases 1 and 2, respectively. The figure shows 
that the release of transferred thermal energy approaches zero, 
which indicates that the transferred thermal energy can be ig- 
nored in the two cases. The irradiation energy plays an important 
role in the early stage of mass overflow and attenuates at the sub- 
sequent stage, which can be explained by the luminosity of the 
primary decreasing rapidly and the luminosity of the secondary 
increasing with mass transfer gradually. The transferred poten- 
tial can exceed the irradiation energy as the mass transfer rate 
grows. The total accretion luminosity in case 2 is higher than the 
one in case 1 because the potentials in the two cases are differ- 
ent. From panel (b), it can be seen that the curve of the accretion 
luminosity fluctuates, indicating that the mass transfer process is 
unstable. 



The total H-burning energy-generation rates of the primary 
in the two cases are shown in Fig. 4. Panel (b) shows the H- 
burning energy-generation rate in case 2 after the main sequence. 
From the difference between curves in panel (a), it is noticed 
that the effect of rotation causes the total H-burning energy- 
generation rate lower. As a result, the evolutive time in the 
main-sequence stage gets longer (cf. Table 1). Moreover, the 
larger fuel supply and lower initial luminosity of the rotating 
stars help to prolong the time which they spend on the main 
sequence (Heger & Langer 2000b). The lifetime extension in 
rotating binary star at the main-sequence stage can also be il- 
lustrated according to Suchkov (2001). Their results show that 
the age-velocity relation (AVR) between F stars in the binary 
system is different from the one between "truly single" F stars. 
The discrepancy between the two AVRs indicates that the pu- 
tative binaries are, on average, older than similar normal sin- 
gle F stars at the same effective temperature and luminosity. 
It is speculated that this peculiarity comes from the impact of 
the interaction of components in a tight pair on stellar evolu- 
tion, which results in the prolonged main-sequence lifetime of 
the primary F star. Moreover, no central helium-burning stage 
exists for case 2 (cf. Table 1). From panel (b), it can be seen 
that the energy-generation rate of the primary vibrates at the 
H-shell burning stage in case 2. These facts suggest that the 
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Fig. 5. Time-dependent variation in luminosity and equivalent radius of primary in two cases. The solid and dotted curves have the 
same meaning as in Fig. 4. 




Fig. 6. Time-dependent variation in surface helium in two cases. The solid and dotted curves have the same meaning as in Fig. 4. 



burning of H-shell is unstable in case 2. The reason lies in the 
centrifugal force reducing the effective gravity at the stellar en- 
velope. The luminosity and surface temperature there decrease 
(Kippenhahn 1977; Langer 1998; Meynet and Maeder 1997). 
Thus, the shell source becomes cooler, thinner, and more de- 
generated as the He core mass increases. As the hydrogen shell 
becomes instable, the thickness y and surface temperature are 

~ 0.203 and 1.1885 x \0 A K, respectively. This physical condi- 
tion leads to thermal instability (Yoon et al., 2004), and the H- 
shell source experiences slight oscillation. It is well known that 
the energy-generation rate is proportional to temperature and 
density (e oc pT")\ therefore, the curve of the H-shell energy- 
generation rate fluctuates. 

The time-dependent variation in the luminosity and the 
equivalent radius of the primary in the two cases are illustrated 
in Fig. 5. Because the rotating star has a lower energy-generation 
rate, the luminosity of the primary is lower, which is the conse- 
quence of decreased central temperature in rotating models due 
to decreased effective gravity (Meynet and Maeder 1997). Then, 
the primary expands slowly in case 2. It is observed that case 1 
reaches point b at t — 2.6725 x 10 7 yr, while case 2 reaches point 
b at / = 2.6267 x I0 7 yr. The initiation time of mass transfer 
for case 2 is advanced by about ~ 1.71%. Similarly, numerical 
calculation by Petrovic et al. (2005b) shows the radius of the ro- 
tating primary increases faster than that of the non-rotating pri- 
mary due to the influence of centrifugal forces. Their results also 



show that mass transfer of Case A starts earlier in rotating binary 
system, which is consistent with ours. If the rotating star is still 
treated as a spherical star, the initiation time of mass overflow 
should be later than that in the non-rotational case. Actually, be- 
cause of the distortion by rotation and tide, the time for mass 
overflow may be extended. Therefore, it is very important to in- 
vestigate distortion in close binary systems. 

The time-dependent variation in the helium compositions 
at the surface of the primary is illustrated in Fig. 6. The H- 
shell burning begins at t = 2.6854 x 10 7 yr in case 1 while at 
t = 2.6329 x 10 7 yr in case 2 (cf. Table 1). Therefore, the ini- 
tiation time of H-shell burning is advanced by 1.71 x I0 5 yr. 
Moreover, the helium composition at the surface of the primary 
is 0.28005 1 at point c, suggesting that the diffusion process pro- 
gresses slowly in a rotating star. Cantiello et al.(2007) also in- 
dicate that rotationally induced mixing before the onset of mass 
transfer is negligible, in contrast to typical O stars evolving sep- 
arately; hence, the alteration of surface compositions depends on 
both initial mass and rotation rates. The sample of the OB-type 
binaries with orbital periods ranging from one to five days by 
Hilditch et al. (2005) shows enhanced N abundance up to 0.4 
dex. Langer et al. (2008) have discovered that for the same bi- 
nary system, but with the initial period of six days instead of 
three days, its mass gainer is accelerated to a rotational veloc- 
ity of nearly 500£:m.s _1 , which produces an extra nitrogen en- 
richment from more than a factor two to about 1 dex in total. 
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Because there is no central helium-burning phase for case 2, the 
diffusion process can be neglected in the interior region of the 
primary after the main sequences. 

4. Conclusions 

The main achievements of this study may be summarised as fol- 
lows. 

(a) The distortion throughout the outer layer of the primary 
is considerable. The detailed theoretical models that investigate 
the outer regions of the two components have deviated some- 
what from the lowest approximation of the Roche model. The 
high-order perturbing potential is required especially in the in- 
vestigation of the evolution of short-period binary system. 

(b) The equilibrium structures of distorted stars are actually 
triaxial ellipsoids. A formula describing rotationally and tidally 
distorted stars is presented. The shape of the ellipsoid is related 
to the mean density of the component and the potentials of cen- 
trifugal and tidal force. 

(c) The radial components of the centrifugal force and the 
tidal force cause the variation in gravitation. The tangent com- 
ponents of the centrifugal force and the tidal force cannot be 
equalized and, instead, they change the shapes of the compo- 
nents from perfect spheres to triaxial ellipsoids. Mass loss and 
associated angular momentum loss are anisotropic in rotating bi- 
nary stars. Ejection is intensified by tidal effect. The ejection of 
an equatorial ring may be favoured by both the opacity effect and 
the higher temperature at the peak of semi-axis b. This effect is 
called the g e (6, ^-effect in this paper. 

(d) The rotating star has an unstable H-burning shell after 
the main sequence. The components expand slowly due to their 
lower luminosity. If the components are still treated as spherical 
stars, some important physical processes can be ignored. 
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